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. Abstract. The application of k t - factorization supplemented with the CCFM small- 

l/"") ' x evolution equation to heavy quark production is discussed. The bb production cross 

sections at the TEVATRON can be consistently described using the k t - factorization 
, formalism together with the unintegrated gluon density obtained within the CCFM 

' evolution approach from a fit to HERA F2 data. Special attention is drawn to the 

comparison with measured visible cross sections. The visible measured cross sections 
at HERA are compared to the hadron level Monte Carlo generator Cascade. 
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1. Introduction 



The calculation of inclusive quantities, like the structure function F 2 (x,Q 2 ) at HERA, 
performed in NLO QCD is in perfect agreement with the measurements. However, 
Catani argues, that the NLO approach, although phenomenologically successful for 
^ ■ F 2 (x, Q 2 ), is not fully satisfactory from a theoretical viewpoint, because "the truncation 
of the splitting functions at a fixed perturbative order is equivalent to assuming that 
the dominant dynamical mechanism leading to scaling violations is the evolution of 
parton cascades with strongly ordered transverse momenta" [ |IJ. As soon as exclusive 
quantities like jet or heavy quark production are investigated, the agreement between 
NLO coefficient functions convoluted with NLO DGLAP [ ||, ||, f|, || parton densities 
and the data is not at all satisfactory: large so-called ^-factors (normalization factors) 
[ H) H 3 are needed to bring the NLO calculations close to the data (K ~ 2 — 4 for 
bottom production at the TEVATRON), indicating that in the calculations a significant 
part of the cross section is still missing. 

At small x the structure function F 2 (x, Q 2 ) is proportional to the sea quark density, 
and the sea-quarks are driven via the DGLAP evolution equations by the gluon density. 
The standard QCD fits determine the parameters of the initial parton distributions 
at a starting scale Qo- With help of the DGLAP evolution equations these parton 
distributions are then evolved to any other scale Q 2 , with the splitting functions still 
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truncated at fixed 0(a s ) (LO) or 0(a 2 ) (NLO). Any physics process in the fixed order 
scheme is then calculated via collinear factorization into the coefficient functions C a (|) 
and collinear (independent of k t ) parton density functions: f a (z,Q 2 ): 

° = vj ^C a (-)f a (z,Q 2 ) (1) 
J z z 

At large energies (small x) the evolution of parton densities proceeds over a large region 
in rapidity Ay ~ log(l/x) and effects of finite transverse momenta of the partons may 
become increasingly important. Cross sections can then be k t - factorized [ |HJ into an 
off-shell (k t dependent) partonic cross section a(^,k t ) and a k t - unintegrated parton 
density function T{z, k t ): 

a = [ -d 2 k t a(-,k t )f(z,k t ) (2) 
j z z 

The unintegrated gluon density T{z, k t ) is described by the BFKL [ [TT| , |13| evolution 



equation in the region of asymptotically large energies (small x). An appropriate 
description valid for both small and large x is given by the CCFM evolution equation [ 



HJ, [15|, [17|, resulting in an unintegrated gluon density A(x, k t , g), which is a function 
also of the additional evolution scale q described below. 

In [ 0] Catani argues that by explicitly carrying out the k t integration in eq. (0) one 



can obtain a form fully consistent with collinear factorization: the coefficient functions 
and also the DGLAP splitting functions leading to f a {z,Q 2 ) are no longer evaluated 
in fixed order perturbation theory but supplemented with the all-order resummation of 
the a s logl/x contribution at small x. 

In this paper heavy quark production at the TEVATRON and at HERA is 
investigated using the k t - factorization approach. The unintegrated gluon density has 
been obtained previously in [ |19| from a CCFM fit to the HERA structure function 
F2(x,Q 2 ). All free parameters are thus fixed and absolute predictions for bottom 
production can be made. As both TEVATRON and HERA data can be described, 
this shows for the first time evidence for the universality of the unintegrated CCFM 
gluon distribution. 

First the basic features of the CCFM evolution equation are recalled and the 
unintegrated gluon density is investigated. Then the calculations for bb production 
at the TEVATRON is presented as well as calculations of the visible cross section for bb 
production at HERA. 



2. The CCFM evolution equation 



A solution of the CCFM evolution equation, which properly describes the inclusive 
structure function F 2 (x, Q 2 ) and also typical small x final state processes at HERA 
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Figure 1. Kinematic variables for multi-gluon emission. The t-channel gluon jour- 
vectors are given by ki and the gluons emitted in the initial state cascade have four- 
vectors pi. The upper angle for any emission is obtained from the quark box, as 
indicated with S. 



has been presented in detail in [ [0|. Figure [I] shows the pattern of QCD initial-state 
radiation in a small-x lepto-production process, together with labels for the kinematics. 
According to the CCFM evolution equation, the emission of partons during the initial 
cascade is only allowed in an angular-ordered region of phase space. The maximum 
allowed angle 5 is defined by the hard scattering quark box, producing the heavy quark 
pair. In terms of Sudakov variables the quark pair momentum is written as: 

Pq + Pq = ^(Pp + ZPe) + Qt (3) 

where p e (p p ) are the incoming electron (proton) momenta, respectively and Q t is the 
transverse momentum of the quark pair in the laboratory frame. Similarly, the momenta 
Pi of the gluons emitted during the initial state cascade are given by (here treated 
massless) : 

rp 1 - 

Pi = Vi{p p + iiP e ) +p ti , & = -3, (4) 

svf 

with Vi — (1 — Zi)Xi-i, Xi = ZiXi-i and s = (p e + p p ) 2 being the squared center of mass 
energy. The variable & is connected to the angle of the emitted gluon with respect to 
the incoming proton and Xi and Vi are the momentum fractions of the exchanged and 
emitted gluons, while Zi is the momentum fraction in the branching (i — 1) — > i and pu 
is the transverse momentum of the emitted gluon i. 

The angular-ordered region is then specified by (Fig. [I]): 

£ < f i < ■ ■ ■ < £n < 2 (5) 
which becomes: 



(6) 
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where the rescaled transverse momenta of the emitted gluons is defined by: 

qi = ^/sYi = : PU (7) 



The CCFM equation for the unintegrated gluon density can be written [ 17, [0 



21f as an integral equation: 

A(x, k t , q) = Ao{x, k t , q) + 

[— I — z q)As(q,zq)P(z,q,k t )A(-,k' t ,q 
J z J irq z \z 

with fcj = + (1 — and g being the upper scale for the last angle of the emission: 
q > z n q n , q n > z n _iqVi-i, qi > Qo- Here q is used as a shorthand notation for the 
2-dimensional vector of the rescaled transverse momentum q = q t = p t j{\ — z). The 
splitting function P(z, q, k t ) and the Sudakov form factor A s (q, zq) are given explicitly 
in[H. 



2.1. The unintegrated gluon density 



In [ [19| the unintegrated gluon density xA(x, kf, q) has been obtained from a fit to the 
structure function F 2 (x 1 Q 2 ) |[ 

In Fig. the CCFM unintegrated gluon density distribution as a function of x and 
k\ is shown and compared to 



2% dxG{Xj /i 



(9) 



with sG(x, /i) being the collinear gluon density of GRV 98 [ in LO and NLO. 
The unintegrated gluon density can be related to the integrated one by: 

r? 

xG{x,n)\ ~ / dk 2 t xA{x,k\,q) (10) 
Jo 

Here the dependence on the scale of the maximum angle q is made explicit: the evolution 
proceeds up to the maximum angle q, which plays the role of the evolution scale in the 
collinear parton densities. This becomes obvious since 

q 2 = x 2 g Es = yxgs = s + Q\ (11) 

The last expression is derived by using Pq + pq — x g p p + yp e + Qt, 2 ~ y/x g and 
s = yXgS — Q 2 . This can be compared to a possible choice of the renormalization and 
factorization scale fi 2 in the collinear approach with /i 2 = Q 2 + 4 • m 2 q and the similarity 
between \i and q becomes obvious. 
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Figure 2. The CCFM kt dependent (unintegrated) gluon density [ [II], 22 
at q = 10 GeV as a function of x for different values of k 2 (upper) and as 
a function of k 2 for different values of x (lower) compared to dxG ^f" - with 
xG(x, fi) being the collinear gluon density of GRV 98 [ [2|| in LO and NLO. 



In Fig. |3] the CCFM gluon density integrated over k t according to eq. (|10|) is 



compared to the gluon densities of GRV 98 [ |23j in LO and NLO. It is interesting 
to note that the CCFM gluon density is flat for x — > at the input scale Q — 1 GeV. 
Even at larger scales the collinear gluon densities rise faster with decreasing x than 
the CCFM gluon density. However, after evolution and convolution with the off-shell 
matrix element the scaling violations of F 2 (x, Q 2 ) and the rise of F 2 towards small x is 
reproduced, as shown in [ [Bl Fig. 4 therein] . A similar trend is observed in the collinear 
fixed order calculations, when going from LO to NLO: at NLO the gluon density is less 
steep at small x, because part of the x dependence is already included in the NLO P qg 
splitting function, as argued in [ (TJ . 

From Fig. |] one can see, that the integrated gluon density from CCFM is larger in 
the medium x range, than the ones from the collinear approach. Due to an additional 

| A Fortran program for the unintegrated gluon density xA(x, fcf , q) can be obtained from [ E2| 
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Figure 3. The CCFM gluon density (solid line) integrated over kt as a function 
of x for different values of Q. For comparison the GRV 98 [ 23] gluon density 
in LO (dashed line) and NLO (dotted line) is also shown. 



1/k 2 suppression in the off-shell matrix elements, the gluon density obviously has to 
be larger to still reproduce the same cross section. In addition only gluon ladders are 
considered in the k t - factorization approach used here, which means that the sea quark 
contribution to the structure function F 2 (x, Q 2 ) comes entirely from boson gluon fusion, 
without any contribution from the intrinsic quark sea. One also should remember, that 
the relation in eq.(^) is only approximately true, since the gluon density itself is not a 
physical observable. 



3. bb production at the TEVATRON 



The cross section for bb production in pp collision at a/s = 1800 GeV is calculated with 
Cascade [ 22|, which is a Monte Carlo implementation of the CCFM approach 
described above. The off-shell matrix element as given in [ |TU| for heavy quarks is used 
with m b = 4.75 GeV. The scale \i used in a s (/z 2 ) is set to /x 2 



m 6+PT ( as m (HI 



with pt being the transverse momentum of the heavy quarks in the pp center-of-mass 
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frame. In Fig. |] the prediction for the cross section for bb production with pseudo- 
rapidity \y b \ < 1 is shown as a function of p™ n and compared to the measurement of 
DO [ RJ. Also shown is the NLO prediction from [ 23, taken from [ In Fig. [5] the 
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Figure 4. Cross section 
for bb production with \y b \ < 
1 as a function of p™ m . 
Shown are the DO [ data 
points, the fixed order NLO 
prediction, and the prediction 
of Cascade. 
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Figure 5. Cross section 
for bb production with \yb\ < 
1 and p% m {b) > 6.5 GeV 
as a function of pJji m (b). 
Shown are the CDF data 
points, the fixed order NLO 
prediction, and the prediction 
of Cascade. 



measured cross section of CDF [ § is shown for 3 values of p™ n (6) with the kinematic 
constraint of \y \, \y b \ < 1 and p™ m (b) > 6.5 GeV together with the prediction from 
Cascade and the NLO calculation from [ 25, taken from [ ||]. In all cases the NLO 
calculation used m& = 4.75 GeV and the factorization and renormalization scales were 



set to 



TTln 



+ p\. Both, DO and CDF measurements are above the NLO 



predictions by a factor of ~ 2. The Cascade predictions are in reasonable agreement 
with the measurements. A similarly good description of the DO and CDF data has 
been obtained in [ |26j using also k t - factorization but supplemented with a BFKL type 
unintegrated gluon density. It is interesting to note, that the CCFM unintegrated gluon 
density has been obtained from inclusive F 2 (x, Q 2 ) at HERA. In this sense, the prediction 
of Cascade is a parameter free prediction of the bb cross section in pp collisions. This 
also shows for the first time evidence for the universality of unintegrated CCFM gluon 
distribution. 

In Fig. ^| the x n and k tn distributions of the gluons entering the hard scattering 
process g+g — > b+b (see Fig. |l|) are shown and the predictions from the k t - factorization 
approach (CASCADE ) are compared to the standard collinear approach (here PYTHIA 
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Figure 6. Comparison of x n and k tn distributions of the gluons entering the process 
g + g — > b + b. Shown are the predictions from Cascade representing k t -factorization 
with the CCFM unintegrated gluon density and also from Pythia representing the 
collinear approach supplemented with initial and final state DGLAP parton showers to 
cover the phase space of p t ordered QCD cascades. 



2~7fl with LO gg — > bb matrix elements supplemented with DGLAP parton showers to 
simulate higher order effects). Whereas the x n distributions agree reasonably well, a 
significant difference is observed in the k tn distribution. However this is not surprising: 
the kt factorization approach includes a large part of the fixed order NLO corrections 
(in collinear factorization). Such corrections are gg —>■ QQg, where the final state gluon 
can have any kinematically allowed transverse momentum, which could be regarded as 
a first step toward a non-p T ordered QCD cascade. 

4. bb production at HERA 

In [ [L9[] the prediction of Cascade for the total bb cross section was compared to the 
extrapolated measurements of the HI [ [8| and ZEUS [ || experiments at HERA. Since 
Cascade generates full hadron level events, a direct comparison with measurements 
can be done, before extrapolating the measurement over the full phase space to the total 
bb cross section. ZEUS [ has measured the dijet cross section which can be attributed 
to bottom production by demanding an electron inside one of the jets. In the kinematic 
range of Q 2 < 1 GeV 2 , 0.2 < y < 0.8, at least two jets with E J T etl{2) > 7(6) GeV and 
\rf et \ < 2.4 and a prompt electron withp^ > 1.6 GeV and \rf \ < 1.1, ZEUS quotes 
the cross section as: 

a b e T^ e++dijet+e - +x = 24.9 ± 6.4+ 4 7 ; 2 pb (ZEUS [ 9]) (12) 
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with the statistical (first) and systematic (second) error given. Within the same 
kinematic region and applying the same jet algorithm Cascade predicts: 

a e+p^e++dijet+e- +X = 20.3^^g pb, (13) 

where the error reflects the variation of m;, = 4.75 =F 0.25 GeV. This value agrees with 
the measurement within the statistical error. The bb cross section, extrapolated from 
the measured jet cross section to the region p\ > 5 GeV and |?7 b | < 2 using Cascade 
is given by: 

a(ep -> e'bbX) = 1.07 ± 0.27 t° f nb (ZEUS using Cascade ), (14) 

which can be compared with the Cascade prediction a(ep — > e'bbX) = 0.87 ± 0.08 nb 
using mb = 4.75 =F 0.25 GeV. The NLO prediction is o~(ep — > e'bbX) = 0.64 nb. When 



using Herwig [ |28j for the extrapolation, ZEUS quotes an extrapolated cross section 
of[|: 

a = 1.6 ± OA(stat.)t° i(syst.)1:oi(ext.) nb 

The quoted extrapolation uncertainty includes the extrapolation with Pythia, which 
is similar to the one obtained with Cascade. Thus the extrapolated cross section 
includes large model uncertainties and the comparison of extrapolated cross sections 
with predictions from NLO calculations are questionable. This also shows the advantage 
of a full hadron level simulation in form of a Monte Carlo program over a fixed order 
parton level prediction. 

It is interesting to note, that the Cascade results agree well with the Pythia 
results for the di-jet plus electron cross section, if heavy quark excitation is included, 
as well as with the extrapolation factor. As in the case of bb production at the 
TEVATRON higher order QCD effects are important, which are already included in 
the k t - factorization approach. 

In Fig. within the same kinematic range the differential cross section for heavy 
quark decays (charm and bottom) as a function of x 7 predicted by Cascade is compared 
with the measurement of ZEUS [ [J. A similar trend as in the case of charm photo- 
production is observed, namely a significant fraction of the cross section with x 7 < 1. In 
LO in the collinear factorization approach this is attributed to resolved photon processes. 
However, in k t - factorization, the x 1 distribution is explained naturally because gluons 
in the initial state need not to be radiated in a pt ordered region and therefore can give 
rise to a high p t jet with transverse momentum larger than that of the heavy quarks. 

The prediction of Cascade has also been compared to the measurement of HI [ 
U for electro-production cross section in Q 2 < 1 GeV 2 , 0.1 < y < 0.8, pj_<2 GeV and 
35° < < 130°: 

a(ep -> e'bbX -> fxX') = 0.176 ± 0.016(siat.)l°;°^(s?/st.) nb 
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Figure 7. The differential cross section do~/dx° bs for heavy quark decays as measured 
by ZEUS [^J and compared to Cascade. Shown are the charm, bottom and the sum of 
both contributions to the cross section (Note there is no additional K factor applied). 



This cross section already includes the extrapolation from measured jets to the muon. 
In the same kinematic range CASCADE predicts: 

a(ep -> e'bbX -> fiX r ) = 0.066tg;{jo7 nb. 

which is a factor ~ 2.6 below the measurement. It is interesting to note, that the ratio 

q(CASCADE) 0.066 
R{m) = a(NLO) = 0^54 = 12 

is similar to the one obtained for the ZEUS extrapolated measurements: -R(ZEUS) = 1.4. 
The measured cross sections of HI and ZEUS cannot be compared directly because 
different kinematic ranges and decay channels were used. To compare both experiments, 
a ratio Rue is defined: 

j~, ^measured /ir\ 

Rmc = (15) 

Using Cascade, the ratios are: 

R MC (m) = a — red = 2.7 ± 0.25181, (16) 

0~MC 

i? MC (ZEUS) = a — red = 1.2 ± 0.321°-^ (17) 

CMC 

where the error comes only from o~ measure( i- This exercise shows that both experiments 
differ by a factor of more than 2 in the published measurements, when compared to 
Cascade. 
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5. Conclusion 

Bottom production at the TEVATRON can be reasonably well described using the 
/^-factorization approach with off-shell matrix elements for the hard scattering process. 
One essential ingredient for the satisfactory description is the unintegrated gluon density, 
which was obtained by a CCFM evolution fitted to structure function data at HERA, 
showing evidence for the universality of the unintegrated gluon density. The comparison 
with the data was performed with the Cascade Monte Carlo event generator, which 
implements A; r factorization together with the CCFM unintegrated gluon density. 

Measurements of bottom production at HERA are also compared to predictions 
from Cascade. The visible dijet plus electron cross section attributed to 6-production 
as measured by ZEUS could be reproduced within the statistical error. It was pointed 
out, that the extrapolation from the measured to the total bb cross section contains 
large model dependencies. If Cascade or Pythia was used for extrapolation, the 
cross section was found to agree with Cascade and even with NLO calculations within 
the combined statistical and systematic error. 

However, the situation is different with the HI measurement: the visible muon cross 
section is already a factor of 2.7 above the prediction from Cascade. Comparing both 
HERA measurements with Monte Carlo predictions of Cascade, it could be shown 
that both experiments, HI and ZEUS, differ in their measurements by a factor of more 
than two. Further measurements, also differential, are desirable to clarify the situation 
at HERA. 

In general the /ct-factorization approach has now proven to be successful even in 
a kinematic region, where typical small- a; effects are expected to be small. It is the 
advantage of that approach that important parts of NLO and even NNLO contributions 
are consistently included due to the off-shell gluons, which enter into the hard scattering 
process. 
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